High light intensity photoassociation in a Bose-Einstein condensate 
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We investigate theoretically the molecular yield in photoassociation of Bose-Einstein condensed sodium 
atoms for light intensities of the order of and above those applied in a recent experiment. Our results show 
that the rate at which ground state molecules may be formed saturates at high light intensities whereas the loss 
rate of condensate atoms does not. This is caused by the opposing roles of the short and long range pair cor- 
relations present near resonance under the influence of the laser and is crucial for the development of efficient 
photoassociation procedures in a condensate. 



PACS numbers: 03.75.Nt, 03.75.Kk, 33.70.-w, 05.30.-d 

The search for efficient ways for photoassociation (PA) 
fT] of Bose-Einstein condensed atoms to condensed diatomic 
molecules has gained much inertia in the recent past. As with 
\ PA in non-degenerate gases (cf. 01 and refs. therein) growing 
theoretical fs, "Tl and experimental 10, IH attention has been 
paid to the investigation of possible limits to the achievable 
PA rate in condensates. Of particular interest is the question 
at which laser intensities / the rate, which is predicted to in- 
crease linearly with / in the low intensity limit, saturates. 

In the experiment reported in Ref. jsl a single color laser 
light pulse has been applied such that the atomic collisions 
occur close to a photoassociative scattering resonance. From 
the measured fraction of atoms remaining after the pulse on- 
resonance loss rates have been deduced for various peak in- 
tensities of the light. These loss rates were found to agree 
well with the two body theory of atoms scattering under the 
influence of a near resonant photoassociative light field 1 3 1 . In 
particular, the experimental results for the time evolution of 
the atom density follow, to a good accuracy, the rate equa- 
tion n{t) = —Kn?{t) if, on resonance, the rate constant 
K is taken to be Kq — (87rfi,/m)|3ao| where oq is the on- 
resonance cornplex scattering length derived from the scatter- 
ing matrix SBl- 

Recent many body calculations |41 have suggested an upper 
limit to the rate of condensate loss due to photoassociation of 
ifpA — {2Trh/m)£, where £ = n^^^^/2Tr is set by the mean 
atomic distance n~^/^. However, it has been pointed out in 
Ref. 1 5 1 that for the maximum atomic density used in the ex- 
periment one has {h/m)£ ~ 3.8 • lO^^^cm'^s^^ which is of 
the order of the PA rate constant Kq reached at the highest ap- 
plied light intensities. The dependence of the measured Kq{I) 
so far gives no evidence for a saturation at higher intensities / 
' 1^1^ and below the unitarity limit \2i- 

In this article we show that while indeed the formation rate 
for ground state molecules saturates for large / the rate of loss 
from the condensate does not, due to the formation of non- 
condensed atom pairs. The theoryapplied in the present work 
is described in detail in Refs. t^|m|9|l and fully takes into ac- 
count the direct loss of condensate atoms into non-condensed, 
motionally excited states, the production and dissociation of 
diatomic molecules during the laser assisted collisions, as well 
as the spontaneous decay of electronically excited atoms. 
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In the experiment fs-l, a A — 591 nm laser couples the 
open channel of electronic \F — l,nip — —1) ground 
state ^^Na atoms, colliding at near threshold energies, to the 
(J = 1,1^ = 135) ro vibrational level of the closed channel 
excimer potential of atom pairs which asymptotically are in 
the electronic states 3^Si/2 + 3^Pi/2 (cf. Fig. 1 in Ref. |51). 
A rectangular intensity profile was chosen, with ramp and fall 
times shorter than 0.5 /is and a length of up to 400 /is. In the 
absence of the laser the collisions of the ground state atoms 
in the { — 1,-1} open channel are characterized by the back- 
ground potential Vbg('')- The laser couples the atoms in this 
channel to a particular rovibrational state \4>^) of the closed 
channel excimer potential Vc\{r) (cf. Fig. 1 in Ref. In 
the absence of Doppler broadening the coupling of the open 
channel to other closed channel states may be neglected. Also 
coupling to other open channels leading, e.g., to inelastic two 
body loss, is sufficiently weak. In the interaction picture and 
the rotating wave approximation (rwa) we may therefore re- 
duce the Hamiltonian of untrapped interacting atom pairs to 

H2B = -^|bg)(bg| + Cl)(i?,-M(^.,cl| 

m 

Here |bg) — |1,— 1)|1,— 1) is the product of hyperfine 
states associated with the { — 1,-1} open channel and |cl) 
is the superposition of products of atomic states associated 
with the closed channel. The closed channel Hamiltonian 
Hc\[r) = —h^V'^lm + Vr\{r) has been replaced by the op- 
erator Ei,Pu, with the (quasi flCrl) projector = |</)^)(0i/| 
onto the resonantly coupled excimer state, Hc\{r)(f)u{r) — 
Eu4>v{i')- The spontaneous decay of to hot atom pairs 
or ground state molecules is taken into account in the imag- 
inary part of — h{uji, — i'f/2) given by the measured 
lifetime r = I/7 = 8.6ns |5]. uj = 27rc/A is the 
laser frequency. The laser induced coupling is described by 
W{x, t) — {h/2)flQ{t) cxp[i(kx — ujt)], with Rabi frequency 
rio(i) whose time dependence takes into account the enve- 
lope of the light field. In Eq. Q the trapping potential may 
be neglected since it does not vary on the scale of the binary 
interactions. 

The microscopic many body approach Q, |^ includes 
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the dynamics of the mean field ^{x,t) = (■!/'g(x))( as 
well as of the density matrix of the non-condensed fraction 
r,,(x,y,i) - (V'](y)V'.(x))t - ^*{y,t)^{^,t) 5,^6,^ and 
the pair correlation function $y (x, y,t) = {ipj{y)ipi{x))t — 
^'(y, i)'I'(x, t)SigSjg of the bosonic field operators for both, 
atoms asymptotically in the ground and excited electronic 
states e {g,e})- All physical quantities of the atom gas 



may be described in terms of correlation functions the dynam- 
ics of which is determined by an infinite hierarchy of coupled 
differential equations. As shown in Refs. Q-Sl this exact hi- 
erarchy can be truncated in a consistent manner Within the 
two-channel description introduced above the mean field of 
the ground state condensate atoms and and the pair functions 
are, to leading order, determined by the coupled equations: 



z;j-*(x,i)=i?iB(x)vI/(x,i) 

( *bg(x,y,t) 
dt I $ci(x,y,t) 



j dy**(y, t) (Hg(|x-y|) [<i>bg(x,y, t) + vl/(x, <)vl/(y, t)\ +hn{x~y, t)<i>ci(x, y, t)) , (2) 



ff'^(x,y,i) 



<fbg(x,y,t)\ / Vbg(|x~y| 



$ci(x,y,t) / V ftf^(x-y,0 



r 



*(x, i)*(y,t). 



(3) 



Here i?g (x) = -W-^I"^ jlm + Vtrap(x) is the Hamiltonian 
of a single atom in the electronic ground state. r2(x, <) = 
fio(i)cos(wt — kx) describes the laser field. Furthermore 
i?^^(x, y,i;) is the general two channel two-body Hamilto- 
nian of interacting atom pairs which in the interaction picture 
and rwa has the form given in Eq. Q. Inelastic loss phe- 
nomena due to e.g. three body recombination will be 
neglected here. 

Using Eq. Q the closed system of Eqs. 0, may be 
solved to obtain a single non-Markovian non-linear Schrodin- 
ger equation for the mean field of the condensed electronic 
ground state atoms 

zn^*(x,i) =i/iB(x)vI/(x,t) 

(■oo o 
~^*{x,t)J^ dT^\x,T) — h{t,T), (4) 

where /i(f,r) = {2TThf{Q,hg\V2B{t)U2Yi{t,T)\{),hg) e{t - 
t). V2B is given by the third term in Eq. Q, U2B{t, t) is the 
time evolution operator corresponding to the Hamiltonian Q, 
6{t — t) is the step function that evaluates to zero except for 
t > T, and |0,bg) is the zero momentum plane wave state of 
the relative motion of two atoms in the open channel. In deriv- 
ing Eq. we have assumed that at the initial time to the gas 
is effectively uncorrected, e.g., a weakly interacting ground 
state Bose-Einstein condensate, fully described by ^'(x, to)- 

At the low relative kinetic energies in a condensate the bi- 
nary scattering properties in the absence of laser light are 
accurately determined by the scattering length ahg_as well 
as by the energy of the last bound state of Vhg II3ll . For 
the ^'^Na atoms in the {—1,-1} channel we use the theo- 
retically determined value of Obg — 54.9aBohr (aBoiu- — 
0.052918 nm) |14]. The energy of the last bound state (vi- 
brational quantum number 1/ = 65) has been measured to be 
Ehg/h = 317.78MHz | lH. Any model of Vbg which prop- 
erly accounts for Obg and i?bg will recover the correct low 
energy scattering properties. To solve Eq. (0} we apply the 
minimal model introduced in Ref. f?]: Vbg = |Xbg)^bg(Xbg|, 
with m^bg/(47rfi^) = — 191.4aBohr andtheconvenientGaus- 



sian function Xhg{r) = cxp(— r^/2(T^^)/(%/27rcrbg)^, fJbg — 
23.97 aBohr- 

A similar ansatz is chosen for the coupling: ^Qq\(P^) — 
|x)C, with e = -(27r)57GHz (a^^^J (kW /cm^))'^/^ and 
x(r) ~ exp(— r^/2o-^)/(\/27ro-)^, cr = 3.55aBohr- Using 
this parametrization the zero momentum limit of the scattering 
amplitude derived from Eq. Q is equal to 47r?i^a(A, I)/m, 
with a given by (cf. fsj]) 



a(A,/) = flbg 



r(/)/2fc 



fc^o A + (5^(/) -^7/2' 



(5) 



where A = Rc — a; is the detuning of the laser fre- 
quency uj from the transition energy. r(/) is the width of 
the laser induced transition, depending on the light intensity 
/ and, for small fc, linearly on the relative momentum hk 
of the colliding atoms |3]. We use the theoretical value of 
d{T{I)/k)/dI = 0.69743 ms-i(kW/cm)-i 6^1) is the 
intensity dependent shift of the resonance which experimen- 
tally has been found to depend linearly on /, dSi, (/) jdl = 
-164MHz(kW/cm2)-i 

In Fig. ^ the time evolution of the condensate fraction 
according to Eq. is shown for a 100 /is light pulse of 
intensity /(<) which is ramped, between t — Q and t — 
0.5 /is, linearly to its peak value /max = 1.2kW cm~^ (dash- 
dotted line), the maximum intensity reported in Ref. jsl. 
The detuning A ~ (27r) 196.8 MHz is chosen such that 
it compensates for 5^{Iyaa.x) in Eq. (|5}- The initial num- 
ber of atoms in the spherical harmonic trap with frequency 
i^ho ^ 198 Hz is A^(0) = 4.0 • 10^. This is compared 
to the time evolution with a constant scattering length and 
loss rate, shown as the dotted line in Fig. where the lo- 
cal density evolves according to the Gross-Pitaevskii equation 
(GPE), as n(x, t) = n(x, 0)[l + i^(A, /„iax)n-(x, 0)t]"\ with 
K{/\J) = (4Vm)|3a(A,/)|. 

In Fig. we also compare the solutions of Eq. and 
of the GPE for /,„ax = lOkW/cm^. The essential differ- 
ence between the dashed and dash-dotted loss curves becomes 
apparent when considering the temporally local decay 'rate' 
K{t) = -Nit)/ Jdxn'^{x,t) shown in Fig. [TJ5. While, after 
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FIG. 1: (a) Time evolution of the remaining fraction of the 
initial atom number A'^(O) =4.0- 10® during a 100 /is pulse, with 
peak intensities of 1.2kW/cm'^ (solid) and lOkW/cm^ (dash- 
dotted line). The fraction N{t)/N{0) is given on a logarithmic 
scale. Shown are the solutions of the non-Markovian non-linear 
Schrodinger equation (NMNLSE) ^4) and, for comparison, the 
evolution according to the Gross-Pitaevskii equation (OPE), 
n(x, t) = n(x, 0)[1 + Kon{x,0){t - frisc)]"^, trisc = 0.5 ^is, with 
Ko = (87r/i/m)|S>o(A, 7)|. (b) Time evolution of the temporally 
local decay 'rate' K{t) = —N{t)/ J dxn^(x,t) correspoding to 
the cases considered in (a). Time is shown on a logarithmic scale. 

the 100 /KS length of the pulse the final densities are compa- 
rable, there is a fundamental difference in the amount of con- 
densate loss at a certain time. This underlines the significance 
of the non-Markovian character of the dynamic equation @ 
and shows that the condensate loss at high intensities may not 
be described by a simple rate equation. 

One sees that the many body theory which takes into ac- 
count all pair correlations created by elastic two body colli- 
sions does not give a saturation of the condensate loss rate 
above ca. / = IkW cmT^. Fig. |2] shows, for different 
-fmax, the fraction of atoms lost from the condensate (filled 
squares) at the final time ffi„ =100/is of the light pulse. 
For large /max the loss approaches 100%. This loss is com- 
pared to the on resonance decrease of the total number of 
atoms iVtot = / (ix[|^'(x)p+r(x, x)] due to the spontaneous 
decay of (filled diamonds), calculated from ftA'tot = 
-7 / dR|/dr0*(r)$ci(R,r,t)P |'l6'l, R and r being the 
cm. and relative coordinates. Since in the 2-channel picture 
considered here ground state molecules can only be formed 
through this decay the upper limit seen in Fig.|2]implies a lim- 
itation of the molecular yield for large /max- 

Which are the physical reasons behind these limits? De- 
pending on the magnitude of A, during the ramps the laser 
causes the closed channel bound state energy to closely ap- 
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FIG. 2: Loss fractions 1 - iVi(*fln)/A'^c(0), i = 'c', 'tot', of the 
condensate (filled squares) and total atom (filled diamonds) loss at 
the end of the laser pulse, ffin = 100 /is, for different laser intensities 
/max- The fraction indicated by the diamonds corresponds to the 
number of ground state molecules formed. The solid line shows the 
result of a solution of the GPE in local density approximation. The 
inset shows the same quantities on a logarithmic scale. 
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FIG. 3: Population {Ml — 1)/Nl of the open channel component 
of the dressed state wave function as a function of the renormalized 
detuning A// from the resonance, for / = 1.2kW/cm^ (dashed line) 
and / = lO.OkW/cm^ (solid line). The real part of the dressed state 
energy vanishes at A = (27r) 160 MHz and A = (27r) 163 MHz, 
respectively. Due to the finite lifetime of these states the populations 
are reduced to ca. 50% and 98%, respectively, as compared to 100% 
in the case that the closed chaimel state is stable, 7 = 0. 

proach or cross the open channel threshold. A rapid ramp of a 
from a(A, 0) — abg to a particular value near resonance rep- 
resents a strong sudden change in the binary interactions. The 
colliding condensate atoms are converted both, into molecules 
and directly to excited correlated pairs which no longer belong 
to the condensate |8, si. fl2il . 

Let us consider the negative energy dressed state wave func- 
tion of the atoms close to the scattering resonance. Within 
the two channel description introduced above the bound state 
spectrum of the open channel potential (cf. Eq. Q) is modi- 
fied by the laser induced coupling to |(/)^). In the pole approx- 
imation, in which only the weakest bound state |(/)^\) of Vbg 
is taken into account the interaction picture Hamiltonian Q 
possesses, depending on / and A, one or two eigenstates of 
negative energy. The normalised dressed states are given by 
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FIG. 4: Open channel component of the radial dressed state density 
for 7 = 1.2kW/cm^ (dashed line) and / = lOkW/cm^ (solid line), 
and a detuning A = (27r) 192.3MHz and A = (27r) 1631.2MHz, 
respectively, which is chosen such that ReiJd is just below thresh- 
old, r is given on a logarithmic scale. The wave functions have 
been computed using the minimal ansatz for Vbg and fiol'/'i/) 
which agrees with a complete coupled channels state at radii large 
compared to the van der Waals length \(m,C<ilf?Y^^ =45aBohr 
f9]. The extension in r scales with the modulus of the scattering 
length: |a(192.3 MHz,1.2kW/cm^)| = 161 asohr, |a(1631MHz, 
lOkW/cm^)! =863aBohr. (Cf. also Ref. [Oil-) 
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Gbg(i;d)VK0. 
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where Ghg{Ed) is the Greens function of the open chan- 
nel Hamiltonian evaluated at the in general complex two 
channel dressed state energy E^- At time t — be 
normalized to one, such that (fi^ is normalized if Md = 
(1 + {ci)AW^G^,g_{Ed)'^Gx,g_{Ed)W\cj)^)Y/^. One key feature 
is now that the open channel component significantly con- 
tributes not only to \<jK3i) away from resonance but also to 
the near threshold state which causes the resonant behaviour 
of a. This is shown in Fig. |3l The second important prop- 
erty is that this component is, near resonance, far ranged, on 
the order of the modulus |a(A,/)| of the scattering length. 
To illustrate this we show in Fig. |4] the radial density of the 
open channel component (j)^{r), Eq. (|6j, for peak laser in- 
tensities /max = 1.2 kW cm^^ and /,nax = 10 kW cm~^, and 
detunings A = {2tt) 192.3 MHz and A = {2tt) 1631.2 MHz, 
respectively, which are chosen such that the real part of Ed is 
just below threshold. (Cf. also Ref. flT].) 

This shows that diatomic correlations may be formed with- 
out the necessity that the atoms have to approach each other 
too closely. At high intensities the closed channel state con- 
tributes less to the dressed state which leads to the saturation 
of the ground state molecule production through spontaneous 
decay above ca. / =1.0kW/cm^. The difference between 
iVtot and Nc as seen in Fig.|2lincludes atom pairs in the unde- 
cayed dressed state as well as in positive energy excited states 



of their relative motion. Since the laser intensity is ramped 
back to zero at the end of the pulse, these atoms are expected 
to emerge as a burst. 

In conclusion we have shown that the atom loss from 
the condensate fraction due to near resonant photoassociative 
coupling to an unstable excimer bound state is not limited ac- 
cording to the mean interatomic separation in the sample. The 
loss is due to the formation of correlations in elastic two body 
collisions which is enhanced by means of the long range na- 
ture of the near resonant dressed state scaling with the modu- 
lus of the scattering length. Nevertheless the possible yield of 
ground state molecules is limited and saturates at high intensi- 
ties. The time dependence of the loss process may in general 
not be described by a simple rate equation, in particular within 
a period after the beginning of the laser pulse which is com- 
parable to the duration of a near resonant elastic collision. 

It is proposed to determine the short time loss for sig- 
nificantly larger light intensities as well as the amount of 
atoms transferred to excited states and bound as ground state 
molecules. 
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